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ABSTRACT 

The direct current (DC) electric field near the reconnection region has been 
proposed as an effective mechanism to accelerate protons and electrons in solar 
flares. A power-law energy spectrum was generally claimed in the simulations of 
electron acceleration by the reconnection electric field. However in most of the 
literature, the electric and magnetic fields were chosen independently. In this pa- 
per, we perform test-particle simulations of electron acceleration in a reconnect- 
ing magnetic field, where both the electric and magnetic fields are adopted from 
numerical simulations of the MHD equations. It is found that the accelerated 
electrons present a truncated power-law energy spectrum with an exponential 
tail at high energies, which is analogous to the case of diffusive shock accelera- 
tion. The influences of reconnection parameters on the spectral feature are also 
investigated, such as the longitudinal and transverse components of the magnetic 
fleld and the size of the current sheet. It is suggested that the DC electric fleld 
alone might not be able to reproduce the observed single or double power-law 
distributions. 

Subject headings: Sun: flares — acceleration of particles — Sun: magnetic flelds 



INTRODUCTION 



Particle acceleration remains a mystery in solar flares, as well as in cosmic rays from the 
extrasolar systems. Observationally, the energetic particles always show a single or double 
power-law spectra l behavior. Three mechanisms have been proposed (see lMiller et al.lll997l : 
AschwandenI |2002| . for reviews), i.e., the direct current (DC ) electric flel d acceleration (e.g., 
Sturrock IQGsh. turbulence (or stochast ic) acceleration (e.g.. Miller 1998 ). and shock acceler- 
ation (e.g.. iBlandford fc Ostrikeiill978l ). Magnetic reconnection, as the effective mechanism 
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for magnetic energy release in solar flares, was demonstrated to be able to prov ide environ- 



ment s favorable for all the above-mentioned acceleration mechanisms to work (jChen et al. 
2007h . 



DC electric field provi des the simple st and most direct means of accelerating particles 
out of a thermal plasma jHolmanlEoSsI ). which eve n allows for a n analytical solution of 
each particle trajectory under certain assumptions (IBulanovl 119801 ) . In the case that the 
electric and magnetic fields induced by the accelerated particles are negligible, test particle 
simulations provide a simple and valid approach to study the particle acceleration in a 



reconnection-associated electric and magnetic fields. For example, ISakail (Il992l ) showed that 
protons and electrons can be promptly accelerated (within 1 s) up to ~70 MeV and ~200 
MeV, respectively. In order to compare the resul t ing en ergy spectrum of the DC-accelerated 
protons with observations, iMori. Sakai. &: Zhaol (119981 ) conducted test particle simulations 
with a hyperbolic magnetic field and a uniform electric field. The accelerated protons present 
a universal power-law spectrum, f{E) ~ E^^, with the spectral index 6 being ~ 2.0 — 2.2. 
Thereafter, a lot of test particle si mulations have been performed for electrons as w ell, 
through either full orbit cal culations (Zharkova &: Gordovskvvll2004J : [Hamilton et al.ll2005l ) or 
guiding center calculations (IWood fc Neukirchll2005l ). Almost all of these simulations claimed 
to have obtained a power-law energy spectrum, with a spectral index around 2, similar to 
that for protons. It is noted, however, that the energy spectra in many of the published 
simulation results deviate significantly from a power-law profile. Moreover, according to 
the thick-target model, the resulting bremsstrahlung hard X-ray (HXR) emissions in these 
simulations would present a power-law energy spectrum with a spectral index around 1. 
However, the RHESSI observations indicated that the HXR emissions in most solar flares 
possess muc h softer spectra, wi th the spectral indices falling in the range between 3 and 7. As 
discussed by IChen et al.l ( 120071 ) , several factors may be attributed to such a big discrepancy. 
For example, the parameters of the reconnecting current sheet were chosen arbitrarily due 
to our poor knowledge of the physical conditions in the localized reconnection region in solar 
flares. Most importantly, the electric and magnetic fields, which should be coupled with each 
other through MHD equations, were often prescribed independently. Doing so often results 
in the electric and magnetic fields that are not compatibl e with the MHD equation s. Some 
improvements were made in recent years. For example, IWood &: NeukirchI (120051 ) derived 
the electric field via the Ohm's law after assuming the distributions o f the magnetic fi e ld, th e 
velocity field, and the resistivity. A further progress was made by [Hamilton et al.l (120051 ) . 
who deduced the solution of the hnearized MHD equations by assuming that the longitudinal 
component of the magnetic field (S^) and the velocity (v) are small perturbations. In the 
real case, both and v can be quite large, therefore, it is necessary to obtain the electric 
and magnetic fields by directly solving the MHD equations. 
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In order to obtain self-consistent electric and magnetic fields, in this paper we first 
perform 2.5-dimensional MHD simulations of the resistive evolution of a current sheet. When 
a steady state is reached, the electric and magnetic fields are then taken out for test particle 
simulations. The paper is organized as follows. In §2, the numerical method is described, 
including the MHD and test particle simulations. The shape of the resulting energy spectrum 
is discussed in §3, and a parameter survey is conducted showing the influences of various 
parameters on the spectral features in §4. A short discussion is presented in §5. 



2. PROBLEM SETUP AND NUMERICAL METHOD 



2.1. Electric and Magnetic Field Configurations 

In order to obtain self-consistent electric and magnetic fields, we numerically solve the 
following 2.5-dimensional (i.e., d/dz = 0), tirne-dependent, compressible res istive MHD 
equations with a multi-step implicit scheme (lHulll989l : IChen. Fang, fc Hul 120001 ): 



dp 
di 



V ■ (pv) = 0, 



dt 



dv 11 

— + (v ■ V)v + -VP - -j X B = 0, 

at P P 

dib 

— - + V ■ W - ri^^l) = 0, 
ot 

+ V ■ VB, + ■^r -B-Vv,-V ■ (r/Vfi^) = 0, 



dT 2h - l)v 

^ + v.VT+(7-l)TV-v-^l— ^j-j-Cg/p = 0, 

dt p/3o 



(1) 
(2) 

(3) 
(4) 
(5) 



where 7 = 5/3 is the ratio of specific heats, rj is the dimensionless resistivity, that is, the 
inverse of the magnetic Reynolds number. The last term on the left-hand side of equation 
Q is the field-aligned heat c onduction, where Q = V ■ [T^/^{B ■ VT)B/B^], and C is the 



dimensionless coefficient (see IChen et al.lll999l . for details). The five independent variables 
are the density (p), velocity (v), magnetic flux function (ip), the longitudinal component of 
the magnetic field (Bz), and temperature (T); note that the magnetic field B is related to 
the magnetic flux function through B = V x {ipGz) + BzSz. When nondimensionalizing the 
MHD equations, the characteristic density and temperature are fixed to be po = 1.67 x 10~^^ 
kg m~^ and Tq = 10^ K, respectively. While, the length scale (Lq) and the plasma beta (/?o) 
are free parameters. The velocity is normalized by the isothermal sound speed vq, and the 
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time is therefore by tq = Lq/vq. Another time scale, the Alfven transit time ta, is related 
to To by Ta = a//5o/2to. The dimensionless resistivity t] is distributed as 

^ _ f ?7ocos(57rx) cos[2(y - 4)7r], if |x| < 0.1, |?/ - 4| < 0.25, 
\ 0, elsewhere, 
where rjo is also a free parameter. 

As the initial conditions, uniform density (p = 1) and temperature (T = 1) are assumed. 
The initial magnetic field is chosen to be a force-free field with a vertical current sheet located 
along the y-axis, which is written as 



{2w/tt) cos(7rx/2w) — w {\x\ < w), 
— \x\ > w), 



(7) 



5, 



3"^ + cos^(7rx/2w) {\x\ < w), 

Bg {\X\ > W), 



where the half-width of the current sheet is w = 0.1, Bg is the longitudinal magnetic field 
(i.e., the guide component) in the background. 

The dimensionless size of the simulation box is —3 < a; < 3 and < y < 8. Because 
of the symmetry, the calculation is made only in the top right quadrant. The top [y = 8) 
and the right-hand {x = 3) sides are treated as open boundaries. Symmetry conditions 
are applied to the left-hand boundary (x = 0) and the bottom {y = 4) of the simulation 



quadr ant. The detailed description of the MHD simulation can be found in IChen et al. 



fll999h . 



As the localized resistivity sets in, the elongated current sheet dissipates and collapses 
into an X-type magnetic configuration. After ~ Tta, the dynamics of the whole simulation 
region becomes steady, keeping all physical quantities almost invariant with time. Figure [1] 
shows the distributions of the temperature {gray scale), the magnetic field [solid lines), and 
the velocity [vector arrows) for the case with Bg = 1, /Sq = 0.01, and rjo = 0.02 at t = 8rA. 
Note that at the steady state, the z-component of the magnetic field near the reconnection 
X-point is around Bg. The small rectangular box in Figure [T] indicates the region where the 
resistivity does not vanish and where test particles are injected. 

The corresponding electric field (E) is then determined by the Ohm's law, i.e., E = 
rjV X B — V X B, where E is nondimensionalized by the characteri stic value of VqBq, Bq 



is the normalization unit of the magnetic field (see I Chen et al.lll999l . for details). Thereby, 
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the electric and magnetic fields are obtained for the ensuing test particle simulations. As an 
example, in the case of (3 = 0.01 and riQ= 0.01, the electric field at the reconnection site is 
~ 600 V m~^, which is super-Dreicer. 



2.2. Test Particle Approach 

Test electrons are uniformly distributed in the resistive region, which is indicated by the 
rectangle in Figure [TJ Initially, these electrons have a Maxwellian velocity distribution with 
the local plasma temperature, which is superimposed on the local plasma bulk velocity. The 
motion of each electron is then calculated by numerically solving the following relativistic 
Lorentz equations 

^(7mov) = g(E + vx B), (9) 



where x and v are the particle position and velocity vectors, respectively; 7 = l/yl — v'^/c^ 
is the Lorentz factor, ttiq and q are the rest mass and the charge of the electron, respectively; 
E and B are the electric and magnetic fields obtained from the MHD simulations as described 
in the previous subsection. It is noted here that collisions are neglected in the momentum 
equation ([9]). 

The 4th-order Runge-Kutta-Fehlberg (RKF45) scheme is used to solve the above equa- 
tions, where the time step (At) is adaptive. In order to get smooth energy spectra, 2 x 10^ 
test particles are simulated. 

In order to confirm our test particle code, we repeated t he simulations of proton ac- 



celeration that were performed by iMori. Sakai. fc Zhad (119981 ). with the same electric and 



m agnetic fields, i.e., E = (0,0 , Ez) and B = {ay, (3x, B^). We chose the same parameters as 



m 



Mori. Sakai. &: Zhad (119981 ) and found that the spectral indices are about 1.9-2.3, almost 



identical to their result, i.e., 2.0-2.2. 



3. ENERGY SPECTRUM: POWER LAW OR EXPONENTIAL? 



The electrons are accelerated under the action of the DC electric field promptly. At 
the same time, they gyro-rotate around the magnetic field lines. After certain rounds of 



- 6 - 



gyro-rotation, the accelerated electrons drift out of the resistive region, and then propagate 
away along magnetic separatrix layers between the reconnection inflow and outflow. Since 
the acceleration and the propagation of the electrons are 2-fold rotational symmetric, the 
particles going upward and downward are collected together to construct their energy spec- 
trum. The left panel of Figure [2] plots the energy spectrum in a log-log scale for the case 
with Bg = 1, Lo=100 m, /3o = 0.01, and rjo = 0.02. To our surprise, the spectrum, which is 
quite smooth except for the high energy tail, does not show a single power-law profile (note 
that a single power-law profile is manifested as a linear line in the log- log scale). Apparently. 



it loo ks like a double power-law profile, as proposed in our recent review paper (jChen et al. 



20071 ). The spectral profile is then fitted with a double power-law distribution, i.e., a harder 
power-law spectrum f{E) ~ E'^'"^ in the energy range of 2-40 keV {red line) plus a softer 
power-law spectrum f{E) ~ E~^-^ in the energy range of 50-120 keV {blue line). It is noted 
that the spectral index at the high energy tail in our simulation is significantly higher than 
those obtained in the published test particle simulations, and is in the typical range of the 
electron spectral indices derived from HXR spectral observations. 

The HXR observat ions from RHESSI d id find double (or broken) power-law spectral pro- 
files in some flares (e.g., Holman et al. 2003h . and th eoretical models have also been proposed 



to explain this feature (e.g., IZhang &: Wangll2004l ). The profile in the left p anel of Figure 



[2] is, h owever, somewhat different from the broken power- law distribution of iHolman et al. 



(120031 ) in that the spectral profile changes gradually between the two power-la w parts in our 
simulation results, while the transition is abrupt in the observation fISuil 120071 ). 



In order to clarify the spectral distribution of the accelerated electrons in our simulations, 
we re-plot the energy spectrum in a linear-log scale in the right panel of Figure [2l It is seen 
that, except for the significantly enhanced lower energy tail, the spectral profile is almost a 
linear line in a wide energy range, which means that the spectrum is close to be exponential 
in most of the energy range. Considering the low energy tail, we fit the spectral profile with 
the combination of an exponential and a power-law functions, i.e., f{E) ~ E~^e~^^^° . With 
the spectral index 6 = 0.35 and the rollover energy £^o=23.1 keV, the fitted line {pink), as 
shown in the right panel of Figure [2|, is consistent with the simulation result in most of the 
energy range, from 8 to over 120 keV. 

Although it was claimed by many authors that DC electric field in magnetic reconnection 
leads to a power-law energy spectrum of the accelerated electrons, quite often the spectrum 
obtained from test particle simulations deviates significantly from a perfect power-law dis- 
tribution. Several authors already noted that the energy spectrum of the DC-accelerated 
electr ons may be characteri zed as either power-law or exponential, depending on the value 
of Bz ([Hamilton et al.ll2005l ) or the trapping time of the electrons inside the acceleration vol- 
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ume (lAnastasiadis. Vlahos. &: Georgouliall997l ). The test particle simulations in this paper, 
with self-consistent electric and magnetic fields obtained from MHD simulations, indicate 
that the electrons accelerated by the DC electric field of magnetic reconnection present a 
spectrum with the form 



f{E) ~ E 



-S^-E/Eo 



which is similar to those obtained in the diffusiv e shock acceleration wh en the power-law 
behavior is truncated by a variety of effects (see lEUison fc Ramatylll985l . for discussions). 
The same sp ectral profile was also o btained in the electron acceleration by random DC 
electric field ( lAnastasiadis et al.l 120021 ) . It is, however, mentioned here that among all our 
simulated scenarios, only in one case with Bg = 1.0, Lo=50 m, /Sq = 0.01, and rjo = 0.005, 
the energy spectrum does show a clear broken power-law shape. 



4. PARAMETER SURVEY 

Several parameters can affect the energy spectral profile of the DC-accelerated electrons 
in our simulations, where the spectral profile is characterized by 6 and Eq as indicated by 
equation fill I) . The free parameters include the longitudinal component of the magnetic 
field (Bg), the length scale (Lq), the resistivity (r/o), and the magnetic field strength, as 
represented by the plasma beta (/3o)- In the following subsections, the effect of varying each 
parameter is investigated individually with other parameters keeping fixed. 



4.1. Effect of the Guide Field B„ 



By theoretical analysis of electron acceleration in a reconnecting current sheet. iLitvinenko 



( 119961 ) propounded that a longitudinal component of magnetic field is necessary to explain 
the accelerated electrons with energy up to 100 keV. In our simulations, the longitudinal 
component (i.e., the guide component along the z-direction) of the magnetic field near the 
reconnection X-point is around Bg. Therefore, we calculated six cases with Bg increasing 
from to 1.0, while the other parameters are Lo=50 m, (3q = 0.01, and r/o = 0.02. The 
energy spectra of the accelerated electrons in the six cases are depicted in the upper panel 
of Figure [31 It is clear that, as Bg increases, more and more electrons are accelerated to 
high energies, which confirr ns the importanc e of longitudinal magnetic field in accelerating 
electrons to higher energies (jLitvinenkdll996l ) . The dependence of S and Eq on Bg is plotted 
in the lower panel of Figure [31 It is found that as Bg increases, Eq increases steadily, while 
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6 first falls rapidly when Bg increases from to 0.4, and then increases slowly. 

In order to study the efficiency of DC acceleration, we calculate the percentage of 
electrons whose energy is higher than 10 keV, and define it as the acceleration rate. We 
find that when Bg is zero, few particles are accelerated and the acceleration rate is 0.0475%. 
However, when Bg reaches 1, the acceleration rate jumps to 11.3%. 



4.2. Effect of the Length Scale Lq 

The size of the current sheet is an important but unknown parameter in the flaring 
process. Therefore, it is worth studying how this parameter would affect the energy spectrum 
of the electrons. As shown in Figured], the initial size of the reconnection sheet is 0.2Lo x 
0.5Lo, where Lq is the length scale of the MHD simulations. With Bg = 1, (3q = 0.005, 
and riQ = 0.02 being fixed, test particle simulations are conducted in four cases with Lq 
being 10, 50, 100, and 200 m, respectively. The upper panel of Figure H] shows the energy 
spectra of the accelerated electrons in the four cases. It is seen that as Lq increases, more 
and more electrons are accelerated to high energies. This is easy to understand since as Lq 
increases, electrons travel a longer distance in the electric field before they migrate out of the 
acceleration region due to gyration. The dependencies of 6 and Eq on Lq are shown in the 
lower panel of Figure |H It is seen that Eq increases almost linearly with Lq, while 6 initially 
increases rapidly with Lq, and then saturates after Lo=100 m. We calculated the fraction of 
the electrons with energies exceeding 10 keV, and found that the acceleration rate increases 
from 4.49% to 26.3% as Lq increases from 10 to 200 m. 



4.3. Effect of the Plasma beta /?o 



In the Petschek model for magnetic reconnection (lPetscheklll964j ). the electric field near 
the reconnection X-point is proportional to B^, where Bq is the magnetic field strength 
in the inflow region. Therefore, the spectral shape of the accelerated electrons should be 
sensitive to the magnetic field. Here, the magnetic field (Bq) is characterized by the plasma 
beta, i.e., Pq = 2^qPqRTq/Bq, With the choice of po and Tq as in §2, Bq is related to /3o 
by Bq ~ 7/y/Po in units of gauss. With Bg = 1, Lq=50 m, and tjq = 0.02, we simulated 
four cases with different (3o. The upper panel of Figure [5] shows the energy spectra of the 
accelerated electrons in the four cases. As expected, it is seen that as Pq decreases, i.e., the 
magnetic field increases, more electrons are accelerated to high energies, and the spectrum 
at the high energy tail tends to become harder. The dependence of 6 and Eq on Pq is plotted 
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in the lower panel of Figure [5l which shows that as /3o decreases, Eq increases significantly, 
while 6 increases when (3q decreases from 0.1 to 0.01 and then decreases. The acceleration 
rate increases significantly from 0.169% to 15.2% as /Sq decreases from 0.1 to 0.005. 



4.4. Effect of the Resistivity rjo 

The resistivity in a reconnecting current sheet is caused by the local microscopic instabil- 
ity of plasma and it is still an unknown parameter in the fiaring process. In our simulations 
the resistivity is characterized by the free parameter rjo. With = I, Po = 0.01, and 
Lq = 50 m, five cases with different tjq are simulated. Note that r/o = 0.01 corresponds 
to an anomalous magnetic diffusivity of 7 x 10^ f2 m, ~ 10^ times the classic value, which 



is quite similar to the numerical results of particle-in-cell simulations (jPetkaki et al.l 12003 



Karlicky fc Bartal 120081 ). Such an anomalous magnetic diffusivity, in addition to the small 
normalization length scale, makes the magnetic Reynolds number, i.e., the inverse of tjq, be 
~ 100. The upper panel of Figure E] shows the corresponding energy spectra of the acceler- 
ated electrons. It is found that the acceleration rate slightly increases from 10.1% to 11.3% 
as rjo increases from 0.005 to 0.02 and then shghtly decreases to 9.70% as rjo increases to 0.08. 
The weak dependence of the acceleration rate is actually due to the population decrease of 
the electrons above ~ 20 keV and the population increase of the electrons between ~ 10 keV 
and ~ 20 keV as r/o increases. The dependence of 6 and Eq on r/o is depicted in the lower 
panel of Figure [61 It is seen that both S and Eq decreases with the increase of r]o. 



5. DISCUSSION 

X-ray and 7-ray observations of solar flares indicated that electrons and protons are 
instantly accelerated in the magnetic reconnection process. Analyses of HXR spectral data 
revealed that the accelerated electrons present a single power-law energy spectrum in most 
flares or a double power-law spectrum in some other events. The co-existence of normal and 
reverse type III radio bursts strongly suggests that nonthermal electrons are acc elerated near 



the r econnection X-point within a compact region less than 2000 km in length (lAschwanden 



2OO2I ). Near the reconnection X-point, both DC electric fleld and the reconnec tion-associated 



turbulence may c ontribute to the acceleration of these nonthermal electrons (jSturrocklll968 



Chen et al.l 120071 ). The DC mechanism, as the simplest model, has been studied by various 
groups with help of test particle simulations. Although the resulting energy spectrum of 
the accelerated electrons is somewhat power law-like, the spectral indices seem to be much 
smaller than the typical values derived from HXR observations. Therefore, more work is 
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required to fill the gap between test-particle simulations and observations. 

Using self-consistent electric and magnetic fields obtained from numerical simulations 
of 2. 5- dimensional nonlinear MHD equations, we investigated the energy spectrum of the 
electrons that are accelerated by the DC electric field in a reconnecting current sheet. Con- 
trary to many of the previous studies, our test particle simulations indicate that the ac- 
celerated electrons present a power-law spectrum truncated by an exponential high energy 
tail, i.e., f(E) ~ E~^e~ ^^^° , which is similar to the case of diffusive shock acceleration 



( lEUison fc Ramatyl Il985l ) and the case of random DC-electric field acc eleration when the 



electron trapping time is long (lAnastasiadis. Vlahos. &: Georgoulislll997l ). A parameter sur- 



vey is conducted to investigate how various physical quantities affect the spectral profile and 
the acceleration rate. As the energy spectra are fitted with the function f{E) ~ E~^e~^^^° , 
it is found that 

(1) Eo, the rollover energy, increases with larger Bg (the guide magnetic field), larger 
Lq (which characterizes the size of the reconnection diffusion region), smaller /?o (the plasma 
beta), and smaller r]Q (the resistivity); 

(2) 6 increases with larger Lq and smaller rjo, and it saturates as Lq is larger than 100 
m or rjo is very small. The dependence of 6 on Bg or /5o is not monotonic. 6 reaches the 
minimum when Bg is ~ 0.4 times the antiparallel component of the reconnecting magnetic 
field, while 6 reaches the maximum when /Sq is around 0.01; 

(3) The acceleration rate, defined here as the percentage of the electrons that are acceler- 
ated above to 10 keV, increases with larger Bg, larger Lq, and smaller Pq. It does not change 
much with rjo, although the percentage of the electrons above 20 keV steadily decreases with 
increasing 779. 

The truncated power-law energy spectrum obtained in this paper is somewhat similar 
to the double power-law distribution, although the transition between the low and high en- 
ergy tails is not so abrupt in most cases. The research in this paper reveals that, even with 
self-consistent electric and magnetic fields, we still cannot reproduce the observed single or 
double power-law energy spectra of the nonthermal electrons in the framework of DC-electric 
field mechanism. This discrepancy between test particle simulations and observations might 
be reconciled if other effects are included. For example, the turbulence in the reconnection 
site would greatly e nhance the coll ision rate between nonthermal electrons and background 



particles or waves (jWu et al.l 120051 ). In other words, DC-electric field in the reconnection 
site alone may not be able to explain the observed spectral features of nonthermal electrons. 
The combination of DC-electric field and the turbulence may finally lead to single or dou- 
ble power-law spectral distributions, which were derived from HXR observations. This is 
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definitely an issue requiring further investigation. 
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Fig. 1. — Distributions of the temperature {gray scale), the magnetic field {solid lines), and 
the velocity {vector arrows) at t = Sta- The rectangle around the reconnection X-point is 
the site where test particles are injected. 
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Fig. 2. — Energy spectrum of the nonthermal electrons for the case with = 1, Lq = 100 
m, Po = 0.01, and rjo = 0.02 plotted in a log-log scale {left, which is fitted with a double 
power law) and a linear- log scale {right, which is fitted with a power law truncated by an 
exponential tail). 
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Fig. 3. — Upper. Energy spectra of the nonthermal electrons in six cases with different guide 
magnetic field Bg, where other parameters are fixed (Lo=50 m, /5o = 0.01, and rjo = 0.02). 
Lower. Variations of 5 and Eq with Bg. 
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Fig. 4. — Upper. Energy spectra of the nonthermal electrons in four cases with different 
length scale Lo, where other parameters are fixed [Bg = 1, /3o = 0.005, and 7]q = 0.02). 
Lower. Variations of 6 and Eq with Lq. 
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Fig. 5. — Upper. Energy spectra of the nonthermal electrons in four cases with different /3oi 
where other parameters are fixed [Bg = 1, Lq = 50 m, and r/o = 0.02). Lower. Variations of 
5 and Eq with [3q. 
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Fig. 6. — Upper. Energy spectra of the nonthermal electrons in four cases with different rjQ, 
where other parameters are fixed {Bg = 1, Lq = 50 m, and /3o = 0.01). Lower. Variations of 
5 and Eq with 779. 



